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Abstract 

We determine the stability conditions for a radially symmetric noncommutative 
scalar soliton at finite noncommutivity parameter 0. We find an intriguing relation- 
ship between the stability and existence conditions for all level-1 solutions, in that 
they all have nearly- vanishing stability eigenvalues at critical 6m 2 . The stability or 
non-stability of the system may then be determined entirely by the c/> 3 coefficient in 
the potential. For higher-level solutions we find an ambiguity in extrapolating solu- 
tions to finite 9 which prevents us from making any general statements. For these 
stability may be determined by comparing the fluctuation eigenvalues to critical 
values which we calculate. 
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1 Introduction 



There has recently been much interest, particularly from string theorists, in noncommu- 
tative geometry J2] . Part of this interest has focused on the construction of noncommu- 
tative scalar solitons, first addressed by Gopakumar et al. [[3] . Explicit constructions were 
carried out for 9 = oo using an isomorphism between the star product and the simple har- 
monic oscillator (SHO) basis. These solutions were then extrapolated to large but finite 
9 using the perturbed equations of motion, and some specific examples of solutions were 
constructed. Classical stability was addressed only as an order-of-magnitude argument: 
a stable solution at large 9 changes mass eigenvalues by 0(1/6*), allowing one to conclude 
that angular fluctuations (in the form of U(oo) transformations) give rise to instability in 
some solutions while still ensuring stability against radial fluctuations. Existence of the 
solitons was then more thoroughly addressed in and , where it was discovered that a 
potential bounded from below produces a critical 9 below which there exists no solution; 
an unbounded potential will always yield a solution for any nonzero 9. This creates a 
need for detailed stability analysis at finite 9. In this paper we report progress in this 
direction. 

The paper is structured as follows. In section 2 we review the construction of noncom- 
mutative solitons. In section 3 we present the tools necessary to determine radial stability 
of a solution. In section 4 we generalize on the stability behavior for all theories. We find 
that at the last point of existence, the critical value (9m 2 ) c , all level- 1 solutions have al- 
most neutral stability in that their eigenvalue is virtually identical to a critical eigenvalue 
which we calculate. The small deviation away from this critical eigenvalue is a product 
entirely of the corrections to the vacuum solution and thus to leading order is a function 
of the 3 coefficient in the potential. For a negative coefficient (those found in typical 
false-vacuum solutions) the solution will be unstable, whereas for a positive coefficient it 
will be stable. For higher-level solutions we find an ambiguity in extrapolating solutions 
to finite 9 which prevents us from making any general statements. For these stability 
may be determined by comparing the fluctuation eigenvalues to critical values which we 
calculate. 

2 Construction of Noncommutative Solitons 

One begins with a field theory of a single scalar in (2 + 1) dimensions with noncom- 
mutativity in the (complex) spatial directions z and z, 




(2.1) 
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where d 2 z = dx dy. Fields in this non-local action are multiplied using the Moyal star 
product 

(A*B)(z,z) = e^ 8 * d -*'- d *' di) A{z,z)B(z',z')\ z=zl . (2.2) 

The most convenient way to obtain solutions to this theory is taking the 9 — > oo limit, 
rescaling z — > z\/6 } z — *> zv9. The kinetic term is then negligible and we can focus upon 
solutions to V '(</>) = 0. As is now an operator this may appear difficult, but it becomes 
greatly simplified by using the procedure outlined in M whereupon radially symmetric 
0(r) can be written in terms of projection operators in the SHO basis, 

oo 

<P(r) = 5>»(n|, (2.3) 

n=0 

oo 

= ^c„0 n (r). (2.4) 

n=0 

so (0*0) (r) = 0(r). It is then simple to see we can write all radially symmetric solutions to 
the equation of motion in the form of Eq. (|2.3|) , with c n satisfying an equation of motion. 
Angular transformations may then be performed by a U (oo) transformation, which is an 
exact symmetry in the 9 = oo limit. The moduli space of all solutions in this limit can 
then be divided into unique levels. For example, a "level 2" solution can be obtained by 
acting with a U(N)/U(N — 2), iV — > oo symmetry on 

= c(0 o + 0i). (2.5) 

Stability determination is easy: if c is a minimum (local or global), the solution is stable. 
If it is a maximum, it is unstable. 

Upon making 9 large but finite, one introduces the kinetic term 

2-JT 

K = — Tr [a, </>][</>, at]. (2.6) 
9 

This will now break the [/(oo) symmetry and change the equation of motion. Stability 
determination requires examining the effect of angular fluctuations and radial fluctuations. 
Level 1 solution stability was considered in the following terms by || . Angular fluctuations 
may be introduced in the form of fluctuations of U which connect 0o to m , m^O. These 
force the solution to decay to 0o- Now consider radial fluctuations for this solution. Since 
we expect these to change by only 0(1/ 9), it is stable or unstable as before. 

There is, however, a need for general stability analysis, valid at any finite 9. Zhou 
H and Durhuus et al. || have analyzed the conditions which allow soliton solutions at 
finite 9. It is the stability determination at finite 9 that we present here. We consider 
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angular fluctuations by using a kind of Bogomolnyi bound on 0(1/6*) kinetic fluctuations 
as explained by Gopakumar and Headrick p], 



5K((f>) 



> 



2tt 

g 2 

2tt. 
? 



(Tr0) + 2]T 



a)eH±,\i)eH 



\{a\a\i)\' 



-Tr 



(2.7) 
(2.8) 



where H\\ is the Hilbert space of projection operators corresponding to {0o, . . . , 4> n -x} 
and 7i± is the Hilbert space of all others. This shows that radially symmetric solutions 
saturate the bound, and it is these solutions we focus upon here. One can then use the 
procedure presented in the next section to determine radial fluctuation stabilty at finite 
9. We should comment on how much we can trust our conclusions, if we can only trust 
angular stability to 0(1/ 9). First, angular fluctuations typically cost more energy and 
should be less significant than radial fluctuations at any 9. Second, since it was shown in 
H that (9m 2 ) c > 13.92 for 4 theories, l/9m 2 is likely still a somewhat small parameter. 



3 Radial Fluctuation Stability 

3.1 Determination of Radial Fluctuation Eigenvalues 

One of the nice features of the SHO basis expansion of the soliton solutions is that we 
have an explicit basis {c n } from which to compute energy E({c n }), equations of motion 
dE/dc n = 0, and fluctuation matrix d 2 E/dc n dc m . The latter must have all positive 
eigenvalues for a stable solution. 



The energy is given as 

oo 

E({c n }) = $^[(2n + l)^ - 2(n + l)c n+1 c n + 9V(c n )} 



n=0 



and thus the stability matrix is 

d 2 E 



\9V"(c n ) + 2(2n + l)]5 nm - 2(n + l)(5 n>m+1 + 5 n , m ^) 



which is more clearly written as 



/ 9V"(c ) +2 -2 

-2 9V"(c x ) + 6 -4 

-4 9V"(c 2 )+ 10 -6 

v -e ••• y 



(3.1) 



(3.2) 



(3.3) 
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We may determine the eigenvalues by putting this in the upper-diagonal form 

/ 0V»(co) +2-2 \ 

ev"(c 2 ) + w A - — a — -6 

V ••■ / 

from which we can read off the eigenvalues 



(3.4) 



A„ = 9V"(c n ) + 2(2n + 1) - Ml. (3.5) 

A n -i 

3.2 Restrictions on Eigenvalues 

It seems we have an obvious method to determine the eigenvalues: armed with the {c n } 
of a solution, compute A and see if it is positive. If it is, compute Ai and see if it is, and 
so forth. But we need to determine the positivity of all A n , not just the first few; how can 
we determine all of them? We do this by assuming A n > 0, and then considering what 
bound is placed on A n that ensures A n+ i is positive. This is done by using the recurrence 
relation 

A (2(n + l)) 2 fi v 

n eV»(c n+1 )+2(2n + 3)-\ n+1 ' 



then setting X n+1 = to get 



A , ( 2 (" + !)) 2 (37) 

An> 9V"(c n+1 ) + 2(2n + 3Y [6J) 



Note that this is more restrictive than simply A n > 0. Now assuming that A n+ i > 0, what 
bound on A n now ensures that \ n+2 is positive? Using the relation between A n +i and \ n+ 2 
and then setting the latter to 0, 

K > < 2 '» +1 » 2 , 2 , +2 „, . (3.8) 

eV»[c, M ) + 2(2n + 3) - w^m^W) 

Similarly, this is again more restrictive than the previous condition. Thus the condition 
placed on A n so that all A after it are positive is the obvious limit, 

A > (2(^+l)) 2 (39) 

n eV»(c n+1 ) + 2(2n + 3) - ey ,, ( Jg^ +5) _... 

= \Z lt (3.10) 
(2(n + l)) 2 



9V"(c n+1 ) + 2(2n + 3) - \^\ 



(3.11) 



This appears to be of little practical use since we need to know all {c n } explicitly. But 
since we are only interested in finite-energy solutions we know all c n — * for large n and 



4 



we can get away with V" (c n ) ~ m 2 + 0(c n ). Thus our strategy is this: compute the first 
few X n by hand using the {c n }, making sure each is positive. When the remaining c n are 
small (c n = at 9 = oo), simply check whether the last A n computed by hand exceeds 
X^ lt plus a small correction due to the small c n . If it does, the solution is stable. Since 
we already have explicit expressions for A n , the first part of this is already done; all that 
remains is to calculate explicit expressions for A^"* for V"(c n ) = m 2 and work out the 
correction terms. This can be done by considering a general potential 

allowing us to approximate V"(c n ) ~ m 2 + 2ac n . If a = 0, we use the next non-vanishing 
term. We then compute X n (8m 2 ) and the perturbation AX^' lt (9m 2 ): 

AAf' = (s^^ + AA-^)^' (3.14) 



_x 

2fn 



cvit \ 

f^TyJ (2^c n+1 -AA^). 



(3.15) 



3.3 Calculation of A£ n *(0ra 



We now use Eq. ( |3.13|) to derive an explicit formula for X c ^ lt {9m 2 ). At first glance 



this appears to be just our stability matrix eigenvalues if all V"(c n ) = m 2 , and thus 
Aq" 4 = 8m 2 + 2, etc. In fact this is not the case. The reason is that this has two 
roots, determined by the method by which we compute the A n . Computing in the order 
A , Ai, . . . , A n _i selects one root, whereas computing X^, . . . , A n+ i selects the other. Since 
our A£™* clearly depends upon the higher A's as just explained, we need to choose this 
latter root. In the appendix we show that 

X^ lt {9m 2 ) = 9m 2 + 2 2 a 2 (3.16) 

where T(a, z) is the incomplete Gamma function. This can be used with the recursion 
relation Eq. (|3.13| ) to easily calculate the next terms. We plot X^iOm 2 ) in Figure 1; the 



qualitative behavior for all levels is similar. 



Two important results are immediately obvious from the solution. First, as 9m 2 — > oo 
we see X c ^ lt — > 0, as we would expect from the interpretation that the eigenvalues decouple 
and must simply be positive "on their own" . Second, the critical eigenvalues aren't real 
for m 2 < 0, so, as expected, there are no stable solutions unless there is a stable true 
vacuum. 
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Figure 1: Critical eigenvalue Aq™* as a function of 9m 2 . 



We then consider the corrections A\^ lt , 



AX 



crit 



-29a 



A 



crit 



2 n+1 



Cn+l + 



\ crit 
A n+l 

2(n + 2) 



(c n +2 + 



(3.17) 



Since we know A„ < 2(n + 1 



AA„ | < \29a(c n+1 + c n+2 + 



26a 



(y\c n+l ) + v'( Cn+2 ) + ■■■) 



(3.18) 
(3.19) 



The convergence of Eq. (|3.19|) was proved by so this correction is well-defined. In fact 
the corrections are very small in most circumstances. 



4 General Stability Analysis 

Naturally one wonders what general statements about a solution's stability are possible 
without checking the eigenvalues numerically. 



4.1 Stability of the Gaussian 

Let us focus on the Gaussian case to begin with; this is the unique level-1 solution at 9 = oo 
which is stable with respect to angular fluctuations. The solution requires (e = l/8m 2 ) 

limco -> c, (4.20) 
limc n -> (4.21) 

uniformly in n for n > 1. Let us quickly review the procedure to construct solutions at 
finite 9, given in 0. The perturbed equation of motion is 

{-td 2 + 1)0 = 60o (4.22) 
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for some scale parameter b. This can be solved in Fourier-space as 

**) = («3) 

Using expressions for (j>(k) in terms of Laguerre polynomials L n (x) we can then calculate 
the solution from 

Cn ( e ) = b [°° e ~ XL n( x ) dx _ (4.24) 



jo 1 + 2ex 

The normalization of L n (x) is such that J Q °° e~ x L n ( y x)L m (x)dx = 5 mn . We determine b by 
matching boundary conditions in Eq. Q4.20Q and Eq. (|4.21|) to get 

V'(bF(e)) = m 2 b(F(e) - 1), F(e) = / dx. (4.25) 

J Q 1 + lex 

As explained in more detail by || and |J, the soliton ceases to exist when there is no 
such b. Just as finding solutions to the equation of motion at 9 = oo meant finding the 
real roots of V, we can think of this as finding the real roots of 

V^ ff (bF(e)) = V'(bF(e)) - m 2 b{F{e) - 1). (4.26) 
At e c = l/{6m 2 ) c the solution must also satisfy dVL^/db = (see Figure 2, for example), 

so 

V"(6 c F(e c )) = m 2 (l--^). (4.27) 

Plugging this into our expression for Ao provides a simple answer for the eigenvalue at the 
last point of existence for the solution, completely independent of potential: 

A o ((0m 2 ) c ) = eV"{co{e c )) + 2 (4.28) 

- (^) c (l-^) + 2 (4.29) 

= Xg*({em 2 ) e ) (4.30) 

where we use the identity proved in the appendix. 



Then our Gaussian stability comes down to only determining the sign of AAq™*. If 
positive, it is slightly unstable; if negative, it is slightly stable. Examining the form of 
AAq™* in Eq. ( |3.15| ) we see this is entirely determined by the sign of the <f) 3 coefficient a/3, 
since it was shown in 00 all c n > if they come from c n = at 9 = oo. 



We provide an example of this using the false- vacuum potential given by ||, 

V&) = 30 2 - 1<P 3 + (4.31) 



7 




Figure 2: At (8m 2 ) c , b will be a local minima of V e ff. 



v 




o 



-1 12 3 



Figure 3: A typical false-vacuum potential. 



This is graphed in Figure 3. Choose the solution <\> = 30o at 9 = oo, so that we begin 
with a stable solution. It is easy to show that for this potential (9m 2 ) c = 46.1488. From 
the procedure outlined this implies 

A ((#m 2 ) c ) = 0.0771, (4.32) 
\ c "\{9m 2 ) c ) = 0.0771, (4.33) 
A\° rit ((9m 2 ) c ) = 0.0355. (4.34) 

In Figure 4 we plot Ao and Ag™* near the critical point, seeing that indeed the solution has 
nearly neutral stability at its critical point. When we then consider the effect of raising 
the critical eigenvalue by A\Q rit ((9m 2 ) c ) = 0.0355, we know it is unstable. 

4.2 Higher Levels 

Let us now attempt to extend this to the level-n solution, 

= c(0 o + --- + 0n-i). (4.35) 

Now <j>(x) will obey 

n-l 

{-ed 2 + l)(f)(x) = b j<Pj( x ) ( 4 - 36 ) 

3=0 
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Figure 4: Aq and Aq near the critical point. 



and we can obtain coefficients from 

71-1 



<: »/v,iO (4-37) 

3=0 

where 

roc e ~x 

F l3 {e) = J -^—^-LiWLjix) dx (4.38) 

which is clearly symmetric under i j. The equations of motion combined with boundary 
conditions are 

-y 71—1 

— V'( Ci ) = -h + h F ir (z = 0, . . . , n - 1) (4.39) 

3=0 

Now rewrite this in the {c;} basis using the substitution 

bi = Fr/ Cj (4.40) 

where Fr. is the matrix inverse of F^ = {Fij}. These equations then need to be decoupled 
to obtain solutions for each Cj. For ^ = 00 this decoupling is trivial since Fy(e = 0) = 6ij, 
producing V'(cj) = as expected. But at finite 9 the equations will include terms such as 
V'(V'(ci)) and be of higher order. This implies there will be an increase in the number of 
solutions. It is then potentially ambiguous how to extrapolate 9 = 00 solutions into finite 
9, and thus it is unclear how one could single out the solution which remains at lowest 9. 
To our knowledge, this solution ambiguity has not been discussed before. 



Let us demonstrate this concept explicitly in a level-2 solution. Eq. ( f4.39| ) becomes 



V'(co) = 6d(Foo-l) + 6iF i, (4.41) 

V\c x ) = boF^ + hiFn-l). (4.42) 

Changing basis allows us to write these as 

V'(c ) = aco + Pc u (4.43) 

V\ Cl ) = (3c + 1Cl . (AAA) 
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Figure 5: V^Acq) at several values of 9m 2 . 



where 



a = -(A-F 11 ) 



8 = ^ 
P A ' 
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^(A-^oo) 



A^FooFn-^. (4.45) 



These are easily decoupled to get 



a 



_(V'( Cl )- 7 c 1 )+/5c 1 , 



2 





(V'{c ) - ac ) + 3cq. 



(4.46) 
(4.47) 



We show in Figure 5 the solutions to Eq. ( [4.471 ) at various values of 9, using the same 
false vacuum potential as we did in the Gaussian case. Multiple solutions are clearly 
evident at finite 9, each of which has a different (9m 2 ) c . 
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5 Conclusion 



We have presented the conditions by which radial stability of a noncommutative scalar 
soliton may be determined. One may explicitly calculate the stability eigenvalues of the 
solution and then compare these to critical eigenvalues which ensure the stability of the 
system. For a level-1 solution, we found that the existence and stability of a system 
are closely intertwined, allowing one to make definite statements about the stability of a 
system at its critical existence point. For level-n solitons we do not believe such a general 
statement can be made, due to the nonuniqueness of solutions at finite 9. 

There are many possible extensions of the analysis presented here. The first would 
be a complete analysis of these extra level-n solutions at finite 9. Another would be a 
comparison with the work by Aganagic et al. M in which unstable solitons in noncom- 
mutative gauge theory were addressed. Also interesting would be to place this in the 
context of the large literature in which noncommutative methods have been applied to 
unstable D-branes and tachyonic condensation. Yet another is that our methods only 
address classical stability and it would be very interesting to include quantum stability. 
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7 Appendix: Equivalence of formulae 

Here we present a proof for the equality of 



x(l - l/F{x)) + 2 = \ c ™{x) 



(7.48) 



where 




(7.49) 



(7.50) 



and 




(2(n + l)) 2 



(7.51) 



x + 2(2n + 3) - K+xi x )' 



11 



We can re-write Eq. (|7.48 ) as 



l 2 1 

X + l ro = — r (7.53) 

To prove this we use an identity found in Perron M §82 formula (26) 

x + a + 8 + i (g±iKg±l) (754) 

x + a + p + 6 x+a+f3+5 __ Jo au 
which holds for/3>0,a>0,a;>0. We note that 

poo 

Mm 8 e' u u p ' x du = lhn(3T(8) 

= 1. (7.56) 
Using this we take the a = 0, (3 — > limit of Eq. ( [7.54] ) to find 



(7.55) 



l 2 _1 

x-\ 
1 



X + 1 Tl — = roc e -„ , (7-57) 

(7.58) 



e x T(0,x) 

We thank D. Goldfeld for providing us with this proof. 
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